Electromagnetic total-<i>f</i> algorithm for gyrokinetic particle-in-cell simulations of boundary plasma in XGC

The simplified δf mixed-variable/pullback electromagnetic simulation algorithm implemented in XGC for core plasma simulations by Cole et al. [Phys. Plasmas 28, 034501 (2021)] has been generalized to a total- f that can include, the first time, boundary diverted magnetic geometry with neutral particle recycling, turbulence, and neoclassical physics. implementation is based on pioneering work Kleiber Mischenko [Kleiber al., Phys. 23, 032501 (2016); Mishchenko Comput. Commun. 238, 194 (2019)]. An demonstration performed DIII-D-like, H-mode plasma, including corresponding comparative electrostatic simulation, which confirms necessary higher fidelity understanding of electron heat transport even at low- β pedestal foot vicinity separatrix.